Measurement of rj meson decays into lepton-antilepton pairs 
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A search for rare lepton decays of the rj meson was performed using the WASA detector at 
CELSIUS. Two candidates for double Dalitz decay r] — > e'^e~e'^e~ events are reported with a 
background of 1.3±0.2 events. This allows to set an upper limit to the branching ratio of 9.7x10""" 



(90% CL). The branching ratio for the decay r; — > e+e~7 is determined to (7.8±0.5stat ±0.8si,at) x 
10"'' in agreement with world average value. An upper limit (90% CL) for the branching ratio for the 
rj e'^e~ decay is 2.7 x 10" 
decays is 3.6 x 10"*. 



and a limit for the sum of the 77 M M M M ^-nd rj tt'^tt m M 
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I. INTRODUCTION 

The 7] decays with lepton pairs are closely related to 
the channels with real photons. A direct consequence of 
Quantum Electrodynamics is that a process with a real 
photon should be accompanied by a process where a vir- 
tual photon converts internally into a lepton-antilepton 
pair (fig.^jb). This fact was first pointed out by Dalitz 
in 1951 [l|. The decays can be related to the corre- 
sponding radiative decays with one or two photons using 
Quantum Electrodynamics and by introducing a function 
of the four momentum transfer squared of the virtual 
photons {ql 2)- 92' '^^) ~ the transion Form Factor 

(FF) (an overview is given e.g. in Q). The fo'^ the 
Dalitz decay is equal to the invariant mass squared of 



the lepton-antilepton pair and . 



> 4m^. 



The FF de- 



scribes the structure of the transition region and it is also 
used for the process 7*7* — > 77 where 2 < (space-like 
virtual photons). 

Experimental information is scarce even for not so rare 
•q meson decays with electron-positron pair(s) as seen 
in table HI where measured and predicted branching ra- 
tios (BR) are summarized. Even the branching ratio 



for the 77 — > e^e 7 decay is known with a rather large 
uncertainty (6.0 ± 0.8) x 10"^ 0. It is worth noting 
that the quoted value was obtained as the average of 
two experimental results with a rather large scale factor 
- 1.4. The recent result from the CLEO Collaboration 
(9.4 ± 0.7) x 10"'^ [l^ is larger by three standard devia- 
tions. None of the 77 decays with double lepton-antilepton 
pairs were observed so far. The decays were studied the- 
oretically already 40 years ago by Jarlskog and Pilkhun 
[§[ assuming a FF equal to one. The effect of the FF 
on the BR is expected to be less than 10% for the decay 
rj — > e+e~e"'"e~ [ll|. For decays with M^M~ pair(s) the 
influence is larger since only large values are probed. 

Decays of neutral pseudoscalar mesons into a lepton- 
antilepton pair, P £~^£~, represent a potentially im- 
portant channel to look for effects of physics beyond the 
Standard Model Q. The dominant mechanism within 
the Standard Model is a second order electromagnetic 
process, additionally suppressed by helicity conservation, 
involving two virtual photons P — > 7*7* shown in fig.[T}:. 
Due to the loop appearing in the diagram the decay is 
sensitive to the values of the FF for any ^ of the pho- 
tons in the loop [l3|. The imaginary part of the decay 
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TABLE L The measured and calculated branching ratios for different 77 decay channels with lepton-antilepton pair(s). The 
data are from 01 if not stated otherwise. The upper limits are for 90% CL. Calculations for single and double Dalitz decays 
are from [1 i, O, El ■ 



a) 




FIG. 1: Diagrams for a) single, b) double Dalitz decays of 
a neutral pseudoscalar meson (tt'', 77 or 77') and c) dominat- 
ing conventional mechanism for decay into a lepton-antilepton 
pair. 



amplitude can be uniquely related to the decay width of 
the 77 — > 77 decay. The experimental value of r(77 — > 77) 
leads to a lower limit (the unitarity bound) of the branch- 
ing ratio: BR(77 — > e+e") > 1.7 x IQ-^ when the real part 
of the decay amplitude is neglected 0, This value 
is much lower than for other decays of tt*^ and 77 into 
lepton-antilepton pairs. This makes the rj — > e~^e~ decay 
rate sensitive to a possible exotic contribution. The best 
experimental upper limit for the BR{rj — > e+e^) comes 
from the CLEO II collaboration 4] and is four orders 
of magnitude higher (table [l|. The decays 7r° — > e+e", 
r] — > fJ-'^ fJ.~ and 77 —^ e+e" are also important in or- 
der to estimate long range contribution to the decay 
Kl The loop diagram of the short-distance 

amplitude is sensitive to the presence of a virtual top 
quark and could be used to improve the knowledge on 
the \Vtd\ element of the CKM matrix 

The real part of the amplitude of the 77 — > e^e^ de- 
cay can be estimated using the measured value of 
BR(77 fJ-^fJ-') [IE [13, El- The assumption that the ra- 
tio between Im and Re parts of the amplitudes for the de- 
cays is the same leads to the prediction BR{ri — > e+e-) k, 



(6 ±0.2) X 10^^. A new, unknown process could increase 
the value. Recently the interest in the decays was revived 
due to the observed excess rate of the tt^ e+e" decay 
[lot with respect to the Standard Model predictions [l^l 
what triggered theoretical speculations that the excess 
might be caused by a neutral vector meson responsible 
for annihilation of a neutral scalar dark matter particle 
[2l| . The consequence could be large (even an oder of 
magnitude) enhancement of the 77 — + e^e^ decay rate. 

The plan of this paper is the following: In part 
II, the experiment is described and the data selection 
is presented. In section II. A the 77 Tr'^vr^vr" decay 
where one of the neutral pions decays via tt" e~^e~^ 
(77 n^n^TT^) is presented. The process is used to verify 
the understanding of the detector response for electrons 
and positrons and to provide normalization for the BR 
of Icptonic 77 decays. This is an extension of the system- 
atical studies from a previous publication that used the 
same data sample [^. In section III. A the Dalitz decay 
77 — > e+e-7 is considered and the BR is determined. In 
sections III.B to III.D the results of the search for the 
r; ^ e+e-e+e", 77 — > 7x+/i-/i+/i-, 77 ^ 7r+7r-/x+/x- and 
77 — > e^e- decays are presented. 



II. THE EXPERIMENT 

The experiment was performed at the CELSIUS stor- 
age rin g in Uppsala, using the WASA detector setup 
(fig. [2]) [22|. Protons with a kinetic energy of 893 MeV 
interacted with frozen droplets of deuterium [23]. The 
77 mesons were produced in the reaction pd -^^He 77 close 
to the 77 production threshold. The detection of '^He ions 
in a zero-degree spectrometer (tagging detector) provided 
a clean 77 trigger independent of decay channel [24| . The 
■^He ions were identified and their energy was measured 
which allowed a clean selection of the pd — >^He 77 reaction 
with a background (mainly due to pd ^^HcTTTr reaction) 
of about 1%. The tagging detector provided a few trig- 
gers per second (at a luminosity of 5 x 10 
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FIG. 2: The WASA detector with zero-degree spectrometer using CELSIUS dipoles. 



yielding on average one recorded 77 event per second. 
During the two weeks of experiment (distributed over a 
period of half a year) nearly 3 x 10^ 77 events were col- 
lected. 

The charged 77 decay products were tracked using a 
cylindrical mini drift chamber (MDC), consisting of 17 
layers of thin-walled (25/im) aluminized mylar tubes and 
built around a beryllium beam pipe of 60 mm diameter 
with wall thickness of only 1.2 mm (3.4x10^'^ radiation 
lengths). Since the target deuterium droplets have a ra- 
dius of 17/im (2xl0~^ radiation lengths) the beam pipe 
is the most important source of photon conversion back- 
ground. For example the fraction of e^e^ pairs from 
7y — !■ 77 with external photon conversion in the beam 
pipe to the Dalitz pairs from 77 — > e+e^7 is about 60%. 
This background could be further suppressed by checking 
the reconstructed position of the vertex of a pair. The 
above features of the WASA detector are crucial for the 
investigation of reactions with e+e" pairs. The MDC is 
placed inside of a superconducting solenoid which pro- 
vides a magnetic field of IT. The MDC is surrounded by 
a barrel of plastic scintillators used mainly to define event 
start time for drift time reconstruction in the MDC. An 
electromagnetic calorimeter consisting of 1012 CsI(Na) 
crystals measures the energies of photons and their im- 
pact points. 

Tracks of electrons and charged pions are reconstructed 
in MDC with an efficiency of about 80% if the transverse 
particle momenta are larger than about 20 MeV/c. One 
should stress that even lepton-antilepton pairs with par- 
allel momenta (and thus minimal value of the invariant 
mass) could be efficiently measured in the MDC. The 
track reconstruction algorithm for the MDC used in the 
present studies was based on a global method of pattern 
recognition in which a constant magnetic field was as- 
sumed. The position resolution of the reconstructed ver- 
tex is about 0.05 cm FWHM in the plane perpendicular 



to the beam and 0.7 cm FWHM along the beam. 

In the offline analysis, events with at least two charged 
particle tracks reconstructed in the MDC were required. 
Events with the tracks originating far from the beam 
target interaction region were rejected. Hit clusters in 
the calorimeter, without associated tracks in the MDC 
and with energy deposit larger than 20 MeV were as- 
sumed to originate from photons. Only events containing 
decay particle candidates with balanced electric charge 
were accepted for further analysis. The results on the 
77 —> 7r+7r~e^e~ decay channel were already presented 
earlier Q. Events with a pair of charged decay prod- 
ucts with opposite electric charges can be attributed 
either to the decay channels with two charged leptons 
or to more frequent channels with two charged pions 
(77 — > 7r+7r~7 and 77 — > TT+Tr^Tr"). 

The following variables are used in the further data 
analysis: 

• The invariant mass of a pair of oppositely charged 
particles (Mee). The electron mass is used in the 
calculations. A clear peak at the lowest value is ex- 
pected for e+e~ pairs from Dalitz decays and from 
conversion of real photons in the detector material. 

• The total invariant mass of all reconstructed decay 
products (for example M{e'^e~j) or M{3tt^)). It 
was required to be consistent with the 77 mass. 

• The missing mass of all decay products {MMrj). It 
should, within errors, be equal to the mass of the 
3 He nucleus (2.808 GeV/c^). 

• The ratio between the momentum measured in the 
MDC and the energy of the shower in the electro- 
magnetic calorimeter associated with the charged 
track (Rp/e)- It permits to distinguish between 
and TT^ when the particles reach the electromag- 
netic calorimeter. 



4 



The opening angle between two reconstructed real 



or virtual photons 



or 6^ 



and the relative 



azimuthal angle between the photons (A0^^» or 
A(/)^,-y*). The angles are given in the laboratory 
frame. 

The separation of electrons from pions relics in the end 
on the kinematics of the reactions studied. Due to the 
large mass difference, the energy momentum conserva- 
tion is violated with the wrong mass assignment. For 
example was the contribution of the background from 
pd ^^HeTT+TT" to the final selection of the rj — > e~^e~ re- 
action found to be negligible. Conversely do neither the 
kinematics or other particle identification methods allow 
us to distinguish pions from muons in the studied chan- 
nels. 



Normalization: 



decay 



In order to normalize the branching ratios of the rj me- 
son decays involving an e~^e~ pair, a monitoring process 
is needed to check the reconstruction efficiency for elec- 
trons and positrons. This is specially important since 
the experiment was split into short time slices distributed 
over a longer time. This data sample was analyzed al- 
ready for a previous paper [gj] on the 77 7r+7r^e^e^ 
decay mode where the r] — > Tr+Tr^Tr*^ decay was used for 
the primary normalization. A cross check was done us- 
ing e+e~7 decays assuming the Particle Data Group 
(PDG) 7] value for this BR. To reduce the systemati- 
cal uncertainty and in order to be able to determine the 
BR(e"'"e~7) we used in this paper both 77 vt+tt^tt'' and 
r] — > TT^TT^TT^ decays for the normalization. 

Dalitz decays of at least one of the three tt^'s from 
the 77 — * TT^Tr'^Tr'^ decays provide an abundant data set of 
events with five photons and an electron-positron pair 
- rj ^ TT^TT^TT^,. The Dalitz decay of the 7r° meson has 
been studied in detail both theoretically and experimen- 
tally. To select a data sample of 77 — > tt'^tt^tt^ events we 
required: 



• at least two tracks from particles with opposite 
charges 

• more than three neutral hit clusters in the calorime- 
ter. 

In fig. [3]the experimental Mee distribution for such events 
is plotted. The peak at low masses is attributed to the 
e^e~ pairs from 77 Tr'^vr^vr*' decays with internal or 
external conversion of one of the photons (solid line in 
the fig. [3]). The maximum at larger masses is due to 
rj decays with a tt+tt^ pair, mainly the rj — > 7r+7r"7r'' 
decay. The relative normalization of the decays differs 
by 15% from what is expected from the branching ratios. 
This difference is attributed to the lower reconstruction 
efficiency for electrons and positrons than for charged 
pions in the MDC. 




[GeV/c^] 



FIG. 3: Invariant mass of the e e~ candidates for events 
with more than three neutral hit clusters. Points - the ex- 
perimental data, solid line - MC simulation for 77 — > tt^'tt^tto 
decay, dotted line - background from decays involving a 
misidentified 7r"'"7r~ pair. 




M(e*ey) [GeV/c^] 



FIG. 4; The e e~7 invariant mass for events with at least 
four neutral hit clusters. Full line - all events, dotted line 
- events with e^e~ invariant mass less than 0.1 GeV . The 
photon giving the M(e^e~7) mass closest to the neutral pion 
mass was selected. 
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The identification of the rj tt^tt^tt^j decay channel is 
confirmed by the reconstruction of the invariant mass of 
the e+e~7 system (M(e+e~7)) where the photon leading 
to the mass value closest to the 7r° mass is selected (fig-H]). 
The M{e~^e~j) distribution is peaked at the 7r° mass 
when the M^e < 0.1 GeV/c^ condition is applied. Fig. [5] 
shows the invariant mass of the three tt*^ 's for the events 
where all pion energies are below 0.2 GeV. For the final 
data sample it was required that the missing mass of the 
system of all decay products is in the range 2.5 GeV/c^ to 
3.0 GeV/c^ and the reconstructed emission angle of the 
r] meson is less than 60°. Assuming that all remaining 
events are due to the decays of rj into three neutral pions, 
the total number of rj mesons Nfj is calculated from the 
formula: 



' (1 - (1 - p)^)A BR{-q ^ ttOttOttO) ^ ' 

where No is the number of the observed events 
(fig. O after background subtraction) and p = 
BR{n° ^ e+e-j) = (1.198 ± 0.032)% 0. The prod- 
uct of the detector acceptance and the reconstruction ef- 
ficiency {A = (13.8 ± 2.0)%) was extracted from a MC 
simulation assuming the Vector Meson Dominance model 
Form Factor for the 7r°. The value of BR{r] tt'^tt'^tt'^) 
is precisely known ~ (32.51 ± 0.28)% [7]. The extracted 




M(37cO) [GeV/c^] 

FIG. 5: The M(37r°) distribution for events with three 
neutral pions reconstructed after all cuts. Points - the ex- 
perimental data, solid line - MC simulation for the sum of 
the 77 — > 7r'^7r''7r^ decay and the background, dashed line - 
T] TT^Tr'^Tr" with external conversion. 

N-q = 256000 ± 18000 agrees with the value from our 
previous paper 6] within two standard deviations. Then, 
for the BR normalization in the present paper we use the 
weighted mean value: N„ = 241000 ± 13000. 



III. RESULTS 

A. Single Dalitz decay ry e^e~^ 

Fig. [5] shows the invariant mass distribution of e^e^7 
candidates (M(e+e^7)) selected by the following condi- 
tions: 

• at least two tracks from particles with opposite 
charges 

• for tracks with matched hit clusters in the calorime- 
ter the condition Rp/E < 1-65 was applied 

• Mee < 0.125 GeV/c2 

• a neutral hit cluster with energy deposit larger than 
180 MeV. 

A clear signal at the M(e"'"e~7) around the 77 meson 
mass is seen. The solid line in Fig. [5] represent MC sim- 
ulation of signal and a sum of all background contribu- 
tions. External conversion of one of the photons from 
rj — > 77 decay comprises the most important background 
as discussed in the section [III The MC underestimates 
the detector resolution in the M(e+e~7). The discrep- 
ancy is caused by not optimal calibration of the calorime- 
ter which is difficult to improve since the data taking was 
distributed over a longer time. We have checked the in- 
fluence of the effect on the extracted values of the BR by 
artificially smearing the MC distributions to match the 
experimental data. 

There are no restrictions on the number of low energy 
neutral hit clusters since due to electron or photon in- 
teraction in the calorimeter an additional hit cluster can 
be created. For about 22% of the events, an additional 
low-energy neutral cluster is reconstructed. The pho- 
ton candidate for the rj — > e^e~^ decay was selected by 
the requirement that the A</)^-y* angle is closest to 180°. 
The signature of the rj — s- e^e~7 decay is an energetic 
photon [E^ > 0.18 GeV). The opening angle 0^^^, is dis- 
tributed between 110° and 150° peaking around 130°. 
Fig. [7| shows 6*^-^* versus M(e"'"e~7). A constraint on 
the angle: 100° < O-y^^, <160° together with a condi- 
tion on the overall missing mass for the decay system 
2.65 GeV/c2 < MM-q <2.90 GeV/c^ cleans the data 
sample significantly. This allows to release the condition 
on Rp/E < 1-65 and this increases the acceptance since 
not all e^e— from 77 — > e"'"e~7 decay reach the calorime- 
ter. Finally 729 events with M(e"'"e~7) between 0.40 
and 0.64 GeV/c^ are identified. The total contribution of 
background (mainly from 77 — > 7r"'"7r~7, 77 — > 7r"'"7r~7r'^ and 
77 77 with one of the photons converting into e^e^ 
pair in the detector material) is estimated to 294 ±15 
events. Fig. [5] shows the M(e"'"e~7) distribution after 
applying all selection cuts mentioned above. 
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M(e*e y) [GeV/c^] 




0.6 0.7 0.8 

M(e*e y) [GeV/c^] 



FIG. 6: The M(e+e"7) distribution for events with Mee < 
0.125 GeV/c^ after particle identification. Points - data, solid 
line is the sum of MC simulations of the signal {rj — > e"^e~7) 
and the background. Dashed line - contribution from 
77 7r'''7r~7r'' and t] ~* 7r^7r~7 decays; dotted line - — > 77 
with one of the photons converted into e"'"e~ pair in the de- 
tector material. 
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FIG. 7: 6-,^, vs M{e'^e 7) before cleaning cuts and without 
particle identification: left - MC simulation for ri e"'"e~7, 
right - experimental data candidates. 

B. The decay rj — > e^e~e^e~ 

In a search for the ij e^e~e~^e~ decay, events with 
exactly two positively and two negatively charged par- 
ticle tracks in the MDC were selected. According to 
the simulations 11% of the reconstructed rj — > e+e^e+e^ 
events should fulfill the following criteria: 

• the relative angle between electron and positron in 
both pairs is smaller than 40° 



FIG. 8: The M(e^e~7) distribution after the final selec- 
tion. Points - experimental data, solid line - MC simulation 
of 77 — > e^e~7, dotted line - MC simulation of 77 —> 77 with 
photon conversion in the detector material. 

• the opening angle between the momenta of the two 
e~^e~ pairs is in the interval 110° to 170° 

• the 77 meson emission angle is smaller than 45° 

• the missing transverse momentum is less than 0.3 
GeV/c. 

In the data only two events passed all selection cuts. The 
event display for one of the two candidates is shown in 

fig.El 

The background is estimated to 1.3 ± 0.2 events and 
originates mainly from single Dalitz decay rj —^ e~^e~"f 
with the photon converting into an e~^e~ pair in the de- 
tector material. 

C. The decays 77 jj.'^ fi^ jj.'^ and 77 ~^ -k'^ fi^ 

The decays 77 —^ fj,^ 11^ and 77 7r+7r^/i+/i^ 
have very similar kinematics. In the analysis we have fo- 
cused on the 77 decay but it is not possible 
to distinguish the two in the present analysis. One starts 
with a similar sample of events as for the 77 — > e+e^e^e" 
decay analysis: four tracks from charged particles with 
charge balance. The events with neutral hit clusters of 
energy larger than 20 MeV or with a track in the For- 
ward Detector (detection angle 2°-17°) are rejected. The 
kinematics is checked assuming muon mass for the four 
charged particles. The opening angle between the mo- 
menta of the two muon pairs is required to be in the 
interval 26°-163°. No candidate event for the discussed 
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FIG. 9: (Color online) Event display for an 
77 — > e^e~ e^e~ candidate event. The shaded area in 
the outermost ring represents the projection of the hit 
calorimeter crystals (the size of the crystals and the radial 
position of the front faces are not to scale). The lines rep- 
resent the reconstructed tracks from the pattern recognition 
program. In addition to layers with straws along the beam, 
the MDC includes twisted layers which cause the spread of 
points for forward/backward going tracks. 



decay channels is left for four candidate muons with in- 
variant mass less than 0.625 GeV/c^ and a missing mass 
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FIG. 10: (Color online) Opening angle between electron 
and positron tracks vs M^e for rj e^e~ event sample se- 
lection: scatter plot - data; shaded area - MC simulation 
of the T) ^ e^e~ decay. A cut corresponding to the ellipse 
shown in the figure, selects 71% of the simulated r) — > e^e~ 
events accepted in the plot. 



D. The decay rj e^e 

Events with two tracks from charged particles of op- 
posite charge are considered. The rj — > e~^e~ decay has 
a distinctive signature in the pd — >'^He7y reaction close to 
threshold: the emitted electron and positron have large 
energies {E > 150 MeV), are co-planar with the beam 
and have a large opening angle (about 130°). In fig. [TUl 



the 



invariant mass is presented as a function of the 



opening angle between the electron and positron for the 
whole data set. The region of the simulated signal after 
reconstruction cuts is also shown. There are no events in 
the region where the majority of the rj — > e^e~ signal is 
expected: e~^e~ opening angle in the interval 120°-160°, 
Mee > 0.49 GeV/c^ and fulfiUing the particle identifica- 
tion criteria 0.5 < Rp/E < 1-65. 



E. Discussion 

The results of the experiment are summarized in ta- 
ble ini The values of the branching ratios are presented 
in table IIIII The confidence limits and intervals for 
decays 1-4 were extracted using Feldman and Cousins 
prescription for small signals with background [25 1. 

The systematical errors were obtained by varying cuts 
applied for selection of the channels and comparison with 
Monte Carlo studies including e.g. different assumptions 



Decay mode 






A 


Events 
background 


Events 
observed 


1. 7; — > e^e~ 


e+ 


e 


(11±1)% 


1.3±0.2 


2 


2. rj ^ fi'^ I-l' 


> 




( 5±1)% 


1.4±0.9 





3. rj ^ ■k'^tt' 


> 




( 5±1)% 


1.4±0.9 





4. rj ^ e^e~ 






(36±3)% 


0.4±0.1 





5. rj —> TT'^ir' 


"e" 




(16±1)% 


7.7±2.0 


24 


6. rj —* e^e" 






(16±2)% 


21.0±2.5 


24 


7. rj ^ e^e~ 


7 




(23±2)% 


294±15 


729 



TABLE II: The detector acceptance A (reconstruction effi- 
ciency included), expected number of background events and 
the number of the observed events after all selection cuts. 



Decay mode 


BR 


BR limit 






90% CL 



1. 


V 






e+ 


e 


(2.7_2 7stat i O.^syst 


xlO"^ 


<9.7xl0" 


-5 


2. 


V 




TT'^n 










<3.6xl0'^ 


-4 


3. 


V 




fi+ji 


"M 








<3.6xl0- 


-4 


4. 


V 




e^e~ 










<2.7xl0" 


-5 


5. 


V 




TT^n 


~e~ 




(4.3^j'gstat i OAsyst 


xlO"" 






6. 


V 




e+e" 










<1.6xl0- 


-4 


7. 


V 






7 




{7.8±0.5stat±0.8syst] 


xlO"^ 







TABLE III: Final results for the branching ratios of lepton 77 
decays. 
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on photon energy reconstruction in the calorimeter and 
on the contribution of interaction with rest gas. The 
main sources of the systematical uncertainty are: 

1. Uncertainty on the total number of the 77 mesons 
in our data sample. That value is dominated by 
a limited number of collected 77 tt^tt^tt^ and 
T] — >• tt+tt"??" decays. Moreover the systematical 
errors of the acceptance and reconstruction effi- 
ciency and BR for these channels were taken into 
account. This contributes with 5.4% to the relative 
uncertainty of the BR for all channels. 

2. Uncertainty of the background contribution. 

3. Acceptance and reconstruction efficiency uncer- 
tainty for a given channel. By using simultaneously 
collected data with similar topologies for normal- 
ization, the contribution of uncertainty of track re- 
construction efficiency is partly canceled. 

For example in the case of 77 ^ e+e~7 decay the sys- 
tematical error is dominated by the uncertainty in the 
acceptance {AA/A = 8%) and it is estimated from the 
discrepancy between the data and the MC. 

The extracted signal for the t] — > e+e~7 decay, 435 ± 
27stat ± ISsyst events, leads to BR(77 e+e~7) = (7.8 ± 
0.5stat±0-Ssyst) X 10"^. This is 20% larger than theoret- 
ical estimates. The result is in between the PDG value 
and the latest CLEO result [l3|. 

The attempt to extract a branching ratio from the ob- 
served two rj —>■ e^e~ e^e~ event candidates leads to the 
value (2.7^2'7stati0.1syst)xl0^^ which is in good agree- 
ment with theoretical estimates. However due to non 
negligible background the value is also consistent with 
zero. If instead one assumes that the events are due to 
background, the upper hmit is 9.7 x 10~^ (90% CL). This 
improves slightly the previous limit from CMD-2 Q . The 
background is mainly due to conversion of the photon 
from the -q e~^e~j decay in the beam tube. It could 
be reduced by checking the position of the reconstructed 



vertex or by selecting events with larger invariant masses 
of the e^e~ pairs. This however decreases the acceptance 
significantly and could not be done in the present study. 

The extracted upper limit for BR(77 e+e~) is 2.7 x 
10^^ (90% CL) and is two times lower than the previous 
one from the CLEO II experiment. 

We also report on the first search for the decays 
T] —> fi'^ fi~ fi'^ fi~ and rj t:^ t:^ fi'^ fi^ . Since the de- 
cays can not be distinguished in the present data 
analysis an upper fimit of 3.6 x 10^"* (90% CL) can 
be given for the sum of the decay branching ra- 
tios. Similarly 77 — > 7r"'"7r~e^e~ and 77 — » e~^e~ ^~ 
decays were not distinguished in the previous analy- 
sis of 7/ — > 7r+7r~e~''e~ decay @. However the branch- 
ing ratio of the r/ e"'"e~/x+/i~ decay is expected to 
be three orders of magnitude lower. Assuming that 
BR(77 — > 7r"'"7r~e+e~) is given as the average of theoreti- 
cal predictions, (3.3±0.3)xl0-'' [8,1,10,12], and taking 
into account the other sources of background reported 
in a limit for the BR(77 e"'"e~/x+/i~) to 1.6x10"'' 
(90% CL) is obtained. For consistency reasons we have 
also reevaluated the result on BR(7/ 7r+7r^e+e~) us- 
ing the Feldman and Cousins approach and the improved 
normalization (table IIII[) . 



Acknowledgments 

We are grateful to the personnel at The Svedberg Lab- 
oratory for their support during the course of the exper- 
iment. The financial support from the Knut and Alice 
Wallenberg Foundation (Sweden), the Swedish Research 
Council and the Goran Gustafsson Foundation (Sweden) 
is acknowledged. This work has also been supported by 
BMBF (Germany) (grants 06HH152, 06TU261), by Rus- 
sian Foundation for Basic Research (grant RFBR 02-02- 
16957) and by the European Community Research In- 
frastructure Activity under FP6, Hadron Physics, RII- 
CT-2004-506078 and HPRTCT-1999-00098. 



[1] R. H. Dalitz, Proc. Phys. Soc. A64, 667 (1951). 

[2] L. G. Landsberg, Phys. Kept. 128, 301 (1985). 

[3] R. R. Akhmetshin et al. (CMD-2), Phys. Lett. B501, 191 

(2001), hep-ex/0012039. 
[4] T. E. Browder et al. (CLEO), Phys. Rev. D56, 5359 

(1997), hep-ex/9706005. 
[5] L. Bergstrom, Zeit. Phys. C14, 129 (1982). 
[6] C. Bargholtz et al. (CELSIUS/WAS A), Phys. Lett. 

B644, 299 (2007), hep-ex/0609007. 
[7] W. M. Yao et al. (Particle Data Group), J. Phys. G33, 

1 (2006). 

[8] C. Jarlskog and H. Pilkuhn, Nucl. Phys. Bl, 264 (1967). 
[9] C. Picciotto and S. Richardson, Phys. Rev. D48, 3395 
(1993). 

[10] A. Faessler, C. Fuchs, and M. I. Krivoruchenko, Phys. 
Rev. C61, 035206 (2000), nucl-th/9904024. 



[11] J. Bijnens and F. Borg (1999), hep-ph/0106130. 

[12] B. Borasoy and R. Nissler, Eur. Phys. J. A33, 95 (2007), 

arXiv:0705.0954 [hep-ph]. 
[13] A. Lopez et al. (CLEO) (2007), arXiv:0707.1601 [hep-ex]. 
[14] S. D. DreU, Nuovo Cimento 11, 693 (1959). 
[15] D. Gomez Dumm and A. Pich, Phys. Rev. Lett. 80, 4633 

(1998), hep-ph/9801298. 
[16] G. Isidori and R. Unterdorfer, JHEP 01, 009 (2004), hep- 

ph/0311084. 

[17] M. J. Savage, M. E. Luke, and M. B. Wise, Phys. Lett. 

B291, 481 (1992), hep-ph/9207233. 
[18] L. Ametller, A. Bramon, and E. Masso, Phys. Rev. D48, 

3388 (1993), hep-ph/9302304. 
[19] E. Abouzaid et al. (KTeV), Phys. Rev. D75, 012004 

(2007), hep-ex/0610072. 
[20] A. E. Dorokhov and M. A. Ivanov, Phys. Rev. D75, 



9 



114007 (2007), arXiv:0704.3498 [hep-ph]. 
[21] Y. Kahn, M. Schmitt, and T. Tait (2007), 

arXiv:0712.0007 [hep-ph]. 
[22] J. Zabierowski et al. (CELSIUS/WAS A) , Phys. Scripta 

T99, 159 (2002). 
[23] C. Ekstrom (CELSIUS/WAS A), Phys. Scripta T99, 169 



(2002). 

[24] C. Bargholtz et al, Nucl. Instrum. Meth. A390, 160 

(1997) . 

[25] G. J. Feldman and R. D. Cousins, Phys. Rev. D57, 3873 

(1998) , physics/9711021. 



